Radial and orbital excitations of static-light mesons 
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We present results for the spectrum of static-light mesons from Nf = 2 lattice QCD. These results 
were obtained using all-to-all light quark propagators on an anisotropic lattice, yielding an improved 
signal resolution when compared to more conventional lattice techniques. With a light quark mass 
close to the strange quark, we have measured the splittings between the ground-state S-wave static- 
light meson and higher excitations. We attempt to identify the quantum numbers of the excited 
states in the context of the reduced spatial symmetries of the lattice. 



I. INTRODUCTION 

The physics of heavy-light hadrons plays a crucial role 
in the determination of CKM matrix elements to test 
the Standard Model and search for new physics. In par- 
ticular, systems with one or more bottom quarks are of 
greatest interest but are also the most theoretically and 
technically challenging. Recently, experiments have be- 
gun to probe the excited state spectrum of B mesons and 
signals for a number of new states have been seen. It is 
still unclear if these new states have the expected masses 
and energies or if they are completely new states, not 
predicted by the standard theory. 

Lattice QCD, in principle, offers the possibility of ab 
initio calculations of the relevant hadronic parameters. 
However, mass-dependent errors which arise in the dis- 
cretisation of the Dirac operator have limited the accu- 
racy of heavy-light simulations. Effective theories, which 
can be used to control or eliminate these errors have 
been developed and widely used in numerical simula- 
tions to predict the masses and decay form factors of 
stable ground states to better than 10% precision 
The static quark approximation is the starting point for 
one such approach. It is based on the observation that 
to a good approximation a heavy-light hadron can be de- 
scribed by an idealised system in which the heavy quark is 
infinitely massive. Heavy-light hadrons exhibit approx- 
imate heavy quark spin and flavour symmetries which 
arise because the Compton wavelength of the heavy 
quark is much smaller than the typical hadronic length 
scale {ruQ ^ ^qcd)- In the static limit these symmetries 
are exact and the hadron consists of light degrees of free- 
dom bound to a static colour source. Systematic correc- 
tions to the static limit can be computed order-by-order 
in AgcD/mg, which forms the basis for heavy quark 
effective theory (HQET). This approach is particularly 
suitable for hadrons containing a bottom quark where 
leading order corrections to the static limit, which in- 
clude chromomagnetic interactions and the heavy quark 
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kinetic energy, are at the ten-percent level. Although 
the static approximation is no longer valid for the much 
lighter charm quark the approach is nevertheless use- 
ful. Relativistic simulations at quark masses below charm 
can be combined with a simulation in the static approx- 
imation to allow interpolation to the charm quark mass 
and in fact this method can also be used to determine b 
physics. 

Given the utility of the static-quark approximation it 
may at first appear surprising that this formalism is so 
under-utilised in lattice simulations. However, simula- 
tions in the static quark approximation have suffered 
from one significant drawback - the signal-to-noise ra- 
tio is very poor. Ultimately, this is due to the fact 
that a static quark propagates only in time. While this 
makes the static quark propagator trivial to evaluate, in a 
conventional lattice simulation where one computes only 
the elements of the light-quark propagator from a single 
space-time site to all other lattice sites, information on 
static quark propagation is only obtained from a single 
spatial lattice site. 

Efforts to overcome this problem have focused on two 
areas. One approach is based on the observation that 
the signal-to-noise ratio of a static-light two-point func- 
tion decays exponentially in time with an exponent which 
receives a large contribution from the static quark self en- 
ergy. The Alpha collaboration realised that by carefully 
choosing the discretisation of the static quark action it is 
possible to reduce the exponent, thereby improving the 
signal resolution at large times [2] ■ 

The second approach, which is used here, is to stochas- 
tically estimate the full Dirac propagator [3, 4,, 5, 6, 7, J, 
i, M, [O, [H, in. El, El m . This allows source and sink 
operators for the static-light correlators to be placed at 
every spatial lattice site, yielding a dramatic increase in 
statistics. In this case, it is crucial to optimise the prop- 
agator estimate to avoid the introduction of unnecessary 
noise. In addition to the vast gain in statistics, the use of 
all-to-all propagators allows for complete freedom in the 
choice of interpolating operators since the measurement 
process is decoupled from the evaluation of the propaga- 
tor. Once the propagator estimate has been computed, 
one can construct any spatially-extended interpolating 
field or apply any level of smearing to the quark fields 
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at both the source and smk, without the need for addi- 
tional fermion matrix inversions. Therefore, an efficient, 
practical all-to-all algorithm [l^ allows us to access not 
only ground state resonances, but also excited hadronic 
states. 

Of course, the two improvements can be combined and 
a simulation incorporating both the improved static ac- 
tion and all-to-all propagators may prove to be even more 
precise. This has not yet been tested. 

This paper is concerned with a precision determina- 
tion of the excited-state spectrum of the Bs meson. The 
spectroscopy of the stable ground-state heavy-light me- 
son is already well-studied on the lattice using a variety of 
approaches including the static approximation [T7| . the 
Fcrmilab [18] approach and NRQCD 19]. However, ac- 
curate simulations of the spectrum of radial and orbital 
excitations, in practice, require an all-to-all propagator 
algorithm and as a consequence are not as mature. 

In this study, we compute the excited-state spectrum of 
static- light mesons in two-flavour QCD. Experimentally, 
little is known of the excited-state spectrum of heavy- 
light hadrons and so predictions from lattice QCD are of 
considerable phenomenological interest. Simulations of 
static-light mesons can also be seen as a natural testing 
ground for all-to-all propagator techniques since they re- 
quire only a single light propagator per configuration. It 
is therefore unsurprising that a number of groups have 
attempted to tackle this particular problem. The static- 
light excited-state spectrum in dynamical QCD has pre- 
viously been studied in refs. [13, Ell- 
in this paper we compute the energies of static-light 
mesons where the light quark mass is close to the strange 
quark. Our results are therefore most relevant for the Bg 
spectrum. In this case, experimentally, the ground state 
(pseudoscalar and vector) energies are well established, 
but there exists no conclusive experimental determina- 
tion of the excited states [12] [U • 

The remainder of the paper is organised as follows. 
Section |ll] describes the group theory of the static-light 
spectrum in the continuum and its analogue on the lat- 
tice. The construction of lattice operators, their analysis 
using a variational approach, the all-to-all propagators 
and lattice actions used are also described. Section IIIII 
discusses the results of the dynamical simulations includ- 
ing finite volume effects and the spectrum of states. In 
Section IIVI we address the issue of multiparticle states. 
Section |V] contains a discussion of the results and in Sec- 
tion |Vl] we draw some conclusions. 



II. THE STATIC-LIGHT SPECTRUM 

The heavy-quark spin symmetry in the static limit 
means that mesons which differ only by the spin of the 
heavy quark are degenerate, ie, there is no hyperfine split- 
ting. In this limit the total angular momentum and par- 
ity of the light degrees of freedom (J^) are conserved 
quantities, and it is conventional to use these quantum 



numbers to label the degenerate hyperfine multiplets [23| . 
So, for example, the pseudoscalar (0^) and vector (1~) 

mesons correspond to = ^ . In the constituent 
quark model description of hadrons these are S-wave (or- 
bital angular momentum L=0) channels. Similarly the 
P-wave channels have Jp ~ ^ , | and the D waves 

are Jp = | 1 1 • Note that we have implicitly assumed 
that the mesons contain a static quark, which has posi- 
tive intrinsic parity. In our measurements we exploit the 
heavy-quark spin symmetry to optimise our signals and 
average over all degenerate channels. To facilitate this 
we construct interpolating operators specifically for the 
light degrees of freedom. We then combine the source 
and sink operators with a simple Wilson line to obtain 
a gauge-invariant two-point function with the required 
quantum numbers [l7| . 

A lattice determination of static-light energies is im- 
portant in its own right but in addition, such an unam- 
biguous determination of the spectrum can also be used 
to indicate which quark models most accurately repro- 
duce the static-light physics. Naively, for a given value 
of L, one expects a natural ordering of energies to pre- 
vail. That is, the multiplet with a larger value of Ji has 
a higher energy. However, some quark models predict 
that at higher orbital and radial excitations this ordering 
should reverse [H, [1^ . Whether this reversal in ordering 
does in fact occur, and at what level it occurs, is still 
a matter of debate. This phenomenon was first investi- 
gated in a dynamical lattice simulation in ref. [26j which 
concluded that no reversal of ordering occurs up to and 
including D-wave excitations. 



A. Spatial rotations on the lattice 

The accurate determination of the excited-state spec- 
trum requires a precise understanding of lattice space- 
time symmetries. The angular momentum quantum 
number J, assigned to states in the continuum, labels 
irreducible representations (irreps) of S'0(3), the group 
of proper rotations which leaves the continuum QCD 
Hamiltonian invariant. On a lattice which is isotropic 
in space, this symmetry group is broken to a 24-element 
subgroup - the octahedral point group O. 

If we base our analysis solely on the continuum sym- 
metry group then we will almost certainly misidentify 
some of the excited-state resonances. We therefore com- 
pute correlation functions of operators which transform 
irreducibly under the lattice symmetry group. These op- 
erators have non-zero overlap with an infinite number 
of states, and the continuum quantum numbers of the 
low-lying states can be deduced by observing (near) de- 
generacies across the lattice irreps. 

The full spatial symmetry group is in fact Oh = 0(^C2, 
where C2 consists of the identity and the space inversion 
operator. The elements of C2 commute with the group 
of proper rotations and the irreps of the full symmetry 
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group follow directly from the irreps of O. Therefore we 
first consider the irreps of this subgroup. 

The octahedral point group has five conjugacy classes 
and accordingly five single- valued irreps. These are la- 
beled Ai, A2, E, Ti and T2, and are of dimension 1,1,2,3 
and 3 respectively. Determining the relationship between 
the angular momentum quantum number and these ir- 
reps is relatively straightforward. Restricting the contin- 
uum irreps, labeled by J, to the elements of O gener- 
ates representations for O which are in general reducible. 
These subduced representations can then be decomposed 
into their constituent irreps. 

The number of times a particular irreducible represen- 
tation of O, labeled by the superscript a, occurs in the 
continuum J irrep is given by the formula 

-'f = i^T.-^xi"\i''- (1) 

Nq is the number of group elements, in this case Nq — 
24. The subscript k labels the conjugacy classes of the 
lattice symmetry group, and Uk is the number of group 
elements in a conjugacy class. is the character of the 

conjugacy class in the lattice irrep and xlf'' is the corre- 
sponding character in the subduced representation. By 
applying this formula systematically for a number of val- 
ues of J, it is possible to deduce the angular momentum 
content of the lattice irreps. 

The preceding analysis applies to single- valued, or 
bosonic, representations only. However, interpolating op- 
erators for the light fermionic degrees of freedom trans- 
form according to the double-valued irreps of O. To find 
these, we consider the double-cover group . This 48- 
element group has eight irreps, but five of these coincide 
with the single-valued irreps of O. The remaining three 
are labeled Gi, G2 and H, and these form the double- 
valued irreps of O. Gi and G2 arc both two-dimensional 
representations, and H has dimension four. The angular- 
momentum content of these irreps can be determined us- 
ing the formula in Eq.([T|), where in this case the sum 
extends over the elements of O^, for which Gi, G2 and 
H are single-valued representations. 

The extension of this discussion to the full spatial sym- 
metry group is clear: Oh has twice the number of irreps 
as O. For example, the irrep Gi of O yields the irreps 
Gig and Giu of Oh- The subscripts g and u denote rep- 
resentations which are even (gerade) and odd (ungerade) 
under spatial inversion, respectively. Table [T] lists the 
even-parity irreps of Oh and the quantum numbers of 
their low-lying constituent states. The table illustrates 
an obvious but important point - states which lie in the 
same continuum irrep, but with difi^erent values of Jz (la- 
beling the rows of the continuum irreps) are in general 
divided between the lattice irreps. For example, states 
with quantum numbers j/'* = corresponding to a 
6-dimensional representation of 0(3), appear in both the 
G2g and Hg irreps. In the continuum, in the absence of 
an external field, these states are degenerate in energy; 



Lattice irrep Dimension 



^Zg ^ 2 ' 2 ■■• 
rr 4 3+ 5+ 7 + 
tzl S. 2 '9. '2 •• 

TABLE L Even-parity irreducible representations of the octa- 
hedral point group with the quantum numbers of the lowest- 
lying constituent states. 

however, this degeneracy is broken by lattice artifacts. 
Therefore, in a numerical study one hopes to determine 
the I"*" energy levels by identifying near-degenerate levels 
in the G2g and Hg irreps which converge in the approach 
to the continuum limit. 

We are now in a position to assign the states of interest 
to their respective lattice irreps. We see that the S wave, 
which has J^" — 5 , lies in the Gi„ irrep. The P 
wave and the | P wave appear in the Gig and Hg ir- 
reps respectively. The D-wave multiplets are labeled | 
and I . Both of these appear in the H^ irrep, and the 

I states also arise in the G2U irrep. The lowest-lying 
states in the remaining lattice irrep , the G2g represen- 
tation, are expected to have the quantum numbers 
corresponding to F-wave excitations. 

B. Lattice operators 

To construct operators which transform accord- 
ing to the irreps of Oh, we first identify sets of 
linearly-independent prototype operators which trans- 
form amongst themselves under Oh- The static quark 
is fixed in space and in this discussion it is convenient 
to identify the position of the static quark with the ori- 
gin of a coordinate system. Suitable prototype operators 
then consist of the spin components of the light quark 
field which can sit at the origin, yielding local operators, 
or may be separated from the static quark by a gauge- 
covariant product of link variables. 

The action of Oh on these prototype operators gener- 
ates a representation of the group, the operators are said 
to transform according to the rows of that representa- 
tion. Interpolating operators which transform according 
to the constituent irreps can be obtained by taking linear 
combinations of the prototype operators. A detailed de- 
scription of how one can construct a complete set of basis 
operators from a prototype set for the more complicated 
case of baryon spectroscopy is given in ref. 27]. Here, we 
simply note that by construction, basis operators for a 
given instance of a lattice irrep form an orthogonal set. 

Before we proceed, it is worth clarifying our terminol- 
ogy. Obviously, by choosing different sets of prototype 
operators it is possible to construct a number of basis sets 
for a given irrep, and it may also happen that a particu- 
lar irrep appears more than once in the decomposition of 
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one of the prototype representations. For the remainder 
of the paper, we wiU say that each disjoint set, transform- 
ing according to a particular irreducible representation, 
corresponds to a single 'instance' of that representation. 

We work in a spin basis where the upper components 
of a Dirac spinor have positive parity, while the lower 
components have negative parity. We can therefore re- 
strict ourselves to prototype operators which use only the 
upper or lower components of the light quark fields. Op- 
erators appearing in a given instance of a representation 
will only contain either the upper or lower quark field 
spin components, and we note that exchanging the up- 
per and lower spin components simply changes the parity 
of the operators. 

First we consider local operators. Under Oh the up- 
per components of the quark field transform according to 
the Gig irrep, while the lower components lie in Gi„. All 
other irreps require spatially-extended operators. The 
simplest spatially-extended prototype operators consist 
of light quark field components separated from the static 
quark by straight-line paths of equal length. For a given 
path length, 6 offsets need to be considered - displace- 
ments along the x, y and z axes in positive and neg- 
ative directions. Combining these offsets with two of 
the four light quark field spin components yields a 12- 
dimensional representation. This representation decom- 
poses into Gig, Giu, Hg and All the S, P and 
D-wave states appear in these irreps. The four irreps 
can be studied using only local operators and operators 
with straight-line offsets. However, identification of the 
D-wave states also requires the use of operators which 
project onto the G2u irrep. The simplest such operators 
can be obtained from a prototype set consisting of spin 
components displaced from the origin along diagonals in 
planes spanned by the x, y and z axes. In terms of link 
variables these diagonals are given by the sum of two 
L-shaped paths as shown in Fig. [TJ We can consider pro- 
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FIG. 1: The planar diagonal path connecting the static quark 
to the light quark field components, which is required to con- 
struct operators for G2u, is given by the sum of two L-shaped 
paths. 

totype operators with more complicated paths to try to 
identify operators which better overlap with the states of 
interest. However, for our purposes, operators involving 
these simplest paths yielded good results. In practice, we 
do not use 'bare' links and quark fields in these interpo- 
lating operators. In order to reduce noise in our measure- 
ments, the link variables are stout-smeared 2^29j. In 
addition, as we describe later, Jacobi smearing jsO^ of the 
light-quark fields plays a central role in the determination 
of the excited-state energies. In both cases the smear- 



ing is implemented such that the smeared quark fields 
and link variables have the same transformation proper- 
ties under Oh as their unsmeared counterparts; therefore, 
the statements made above about the choices of path re- 
quired to access each irrep are unchanged. 

Within this collection of operators, there is still some 
freedom in the choice of interpolating operators. This is 
particularly evident in the case of straight-line displaced 
operators, where the prototype representations using ei- 
ther the upper or lower components have the same de- 
composition in terms of the lattice irreps. From the op- 
erator set, final measurements are performed with the 
operators which couple optimally to the states of inter- 
est. The correlation functions we wish to evaluate give 
the amplitude for creating a static-light meson, contain- 
ing a static quark, with the required quantum numbers 
at some initial time and annihilating the meson at a later 
time. When we evaluate the amplitude for this process 
we find that interpolating operators which use the upper 
components of the light quark field yield noisy measure- 
ments. Fortunately, the operators which use the lower 
quark field components have a larger overlap with the 
states of interest and give much cleaner signals. To avoid 
the introduction of unnecessary noise, we restrict our- 
selves to operators containing the lower quark field com- 
ponents in the evaluation of two-point functions which 
describe a meson which contains a static quark, propa- 
gating forwards in time. 

Therefore, the only local operators used transform ac- 
cording to the Giu irrep, and we use spatially-extended 
operators to project onto the Gig representation. As dis- 
cussed, operators with straight-line paths project onto 
the Hg and Hu irreps, and it is possible to construct oper- 
ators incorporating planar-diagonal paths which project 
onto the G2u irrep. Similar operators can be used to 
project onto the G2g irrep. However, the lowest-lying 

states in this representation are the I"*" and I"*" F waves. 
These very heavy states are not of immediate interest to 
us, and we do not include the G2g irrep in this study. 

Therefore it is possible to evaluate correlation func- 
tions for five lattice irreps using only the lower quark field 
components and the simple gauge-covariant paths listed 
previously. However, using the transformation property 
of the fermionic degrees of freedom under time reversal, 
the upper spin components of the light quark field can be 
used to further enhance our statistics. So far, the corre- 
lation functions have been described in terms of a meson 
containing a static quark, propagating forwards in time. 
Formally, the same amplitude can be obtained by revers- 
ing the direction of the Wilson line and substituting the 
upper quark field spin components into the interpolating 
operators. This is a correlation function for the light de- 
grees of freedom with the same angular momentum, but 
opposite parity, propagating backwards in time. It can 
also be regarded as a correlation function for a static- 
light meson containing a static anti-quark. By averaging 
over correlators for mesons containing a static quark, and 
mesons containing a static anti-quark, we ensure that all 
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four light quark spin components are used in our mea- 
surements. 



C. The variational method 

The operators described in the previous section cou- 
ple to all the states appearing in their respective irreps. 
At sufficiently large times it is expected that the be- 
haviour of the corresponding two-point functions will re- 
duce to the usual single exponential decay with an expo- 
nent equal to the energy of the lightest state appearing 
in that irrep. A variational approach must therefore be 
applied to compute the energies of higher-lying states in 
each irrep. To apply the variational method in the lat- 
tice irrep R, we first compute a matrix of correlation 
functions 

C^p{t) cx Y.{n\oi^\t + T,x)W{t + T,t;x)Oj,''\t,x)\n) 

+T.R (2) 

where the subscripts a, (3 label interpolating operators for 
the light degrees of freedom in the irrep R. The barred 
operator contains spin components of the fermion field 
q{x), which lies in the anti- fundamental representation 
of colour SU{3), and the unbarred operator contains the 
spin components of q{x). yV(t -I- T;t,x) is a Wilson line 
from t to t + T at the spatial site x. In Eq. ^ only the 
contribution to the correlator involving the lower quark 
field components is given explicitly, and the initials T.R. 
denote the piece involving the upper quark field compo- 
nents, obtained by time reversal. In terms of the interpo- 
lating operators appearing in Eq. ([2|) , the operator which 
optimally couples to the n*^ excited state is given by 

0i«)(t,x)-^«„„Oi«)(t,a;). (3) 

a 

The real-valued coefficients Vna form the elements of a 
vector v„ which is obtained by solving the generalised 
eigenvalue equation 

C(tl5)v„ = Xn{tD,to)C{to)Vn, (4) 

where to is a fixed initial time and Id is a later 
reference time. The eigenvalues in Eq.Q sat- 
isfy limj^o, A„(t, to) = e-^"(*^*o) [l-KO(e"*"(*-*«))], 
where En is the energy of the n*^ excited state and Sn is 
the absolute value of the energy difference to the nearest 
state. 

The excited state energies can be determined by fitting 
to optimised correlation matrices 

Cn,n{r) CX Y.{n\^l''\t + T,x)W{t + T,t;x^(^\t,x)\n) 
+T.R., 

= v„C(t)v,„. (5) 



To construct the initial correlation matrices, C(t), we 
require a set of operators which transform according to 
the irrep R. However, these operators cannot lie in a 
single instance of R because the correlation between dif- 
ferent rows in a single instance of an irrep is zero, render- 
ing the variational method invalid. To obtain correlation 
matrices we apply a number of levels of Jacobi smearing 
to the quark field components used in the interpolating 
operators which generates several instances of the irrep. 
For each row in R, cross correlations of operators ap- 
pearing in different instances of the irrep are computed. 
This results in a correlation matrix for each row in the 
irrep. These matrices are identical up to statistical errors 
and we average over the rows of the irrep to obtain a fi- 
nal correlation matrix. The variational method (Eq.([4]), 
Eq-®) is then applied to this averaged correlation ma- 
trix to determine the energy of the ground state as well 
as higher-lying states in R. 

D. All-to-all propagators 

The use of an efficient estimate for all elements of the 
lattice Dirac propagator lies at heart of this study. Al- 
though the move from point-to-all fermion propagators 
to all-to-all propagators undoubtedly introduces an ad- 
ditional computational overhead into lattice simulations, 
there exist very many applications in lattice field the- 
ory where the benefit of using all-to-all propagators far 
outweighs this additional cost. The work described here 
involves two areas where all-to-all techniques are particu- 
larly useful. Firstly, the reduction in noise per configura- 
tion achieved by averaging the correlation matrices over 
the lattice volume is essential in studies of static-light 
systems on rather expensive dynamical background con- 
figurations. Secondly, all-to-all propagators are of great 
advantage in excited state spectroscopy. Because all ele- 
ments of the lattice propagator are known, the construc- 
tion of spatially-extended interpolating operators and the 
smearing of the quark fields require no additional inver- 
sions of the fermion matrix. The utility of all-to-all prop- 
agators in constructing operators which optimally couple 
to orbital and radial excitations cannot be overstated. 

To estimate all elements of the Dirac pro pag ator we 
employ the hybrid method introduced in ref. [l6[ . In this 
approach, the contribution of the low-lying eigenmodes 
of the Dirac operator to the fermion propagator is evalu- 
ated exactly, since these modes are thought to dominate 
the long-range physics of interest; the contribution of the 
higher eigenmodes is then estimated stochastically. In 
terms of the hermitian fermion matrix Q — 75^/, the 
propagator is simply 

= Q-I75. (6) 

To proceed, Q^^ is written as a spectral sum which is 
divided into two pieces 

Q-' = Qo + Qi, (7) 
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with 

where Qv^^^ = XiV^"^^ and TV is the rank of the fermion 
matrix. Therefore, according to our method, Qq, which 
contains the first N^^ eigenvectors, is computed exactly. 
Qi is the contribution of the higher-lying modes to the 
propagator, which is estimated stochastically. Using the 
projection operator. Pi — 1—J2i^=i v^^''^v^^^\ this contri- 
bution can be written in a more useful form Qi = Q~^Pi. 

In a naive stochastic estimate, one generates a set of 
Nr random noise sources {r][iy...r]i^pf^^} with the structure 
of quark fields and the property that 

{{V[r]{x)(E>r][r]{yy)) ~ S^y, (9) 

where a;, y denote any set of quark field indices and ((...)) 
denotes an average over the set of noise vectors. An es- 
timate for Qi is then given by 

Qiix,y)~ {{Ar]ix)®vl]))^ (10) 

where the set {"fAfr]} is obtained by solving ■(/'[^J = 
Q^^Piri[r] for each source vector. Unfortunately, the 
variance in this estimate can be quite large, and this error 
decreases as ly/TT^ so that Eq. pH)) becomes exact only in 
the limit N.^ — > oo. However, the stochastic estimate can 
be improved upon substantially through a careful parti- 
tioning or 'dilution' of the noise sources. The idea here 
is to break each original noise source into disjoint pieces 
according to some set of quark field indices. For exam- 
ple, in time dilution, which appears to play a particularly 
important role in variance reduction, each noise vector is 
partitioned as 

Nt-l 

V[r]ix,t) = E r]P^{x,t), (11) 

i=0 

where rj^^^ (x, i) = for i 7^ t. 

One can dilute in any combination of quark field in- 
dices, generating Nd diluted noise vectors from each orig- 
inal source. Inverting the fermion matrix on these diluted 
sources yields an unbiased estimate for from a single 
original noise source 

^lp,t)®rjl^^\xo,to). (12) 

i=0 

In this study, we use noise and each component of 
the noise sources has modulus 1. For such noise, the 
limit of full dilution corresponds to an exact evaluation 
of Qi- This is in stark contrast to the naive estimate 
where an infinite number of matrix inversions is required 
to evaluate the propagator exactly. 

Averaging each diluted estimate over the distribution 
of noise sources means that a total of Ninv = x N,. 



matrix inversions are used in the estimator of Qi. There 
is considerable freedom in how we partition the noise 
vectors and the optimal dilution path is that dilution 
which yields satisfactory results while minimising N^v ■ 

Although the number of eigenvectors, the number of 
source vectors and the dilution level required depend very 
sensitively on the details of the measurement, in our stud- 
ies of static-light systems we have found that computing 
the contribution of a small number of the low- lying eigen- 
modes exactly, and using just one or two noise sources 
with a moderate level of dilution yields excellent results. 
The precise details of the fermion propagator estimates 
used in this study are given in Section IIIFI 

E. Lattice actions 

The lattice gauge and Dirac actions used in this study 
are described in detail in refs. [3l|, [s^l- They are formu- 
lated specifically for anisotropic lattices. The gauge ac- 
tion is Symanzik and tadpole improved with leading dis- 
cretisation errors of ©(a^, , a^a^), where Og is the lat- 
tice spacing in the spatial directions and at is the lattice 
spacing in the temporal direction. The fermion action 
is Wilson-like in the temporal direction but the lattice 
Lagrangian includes an operator of dimension seven to 
lift the spatial doublers. This action has leading discreti- 
sation errors of ©(a^, Ot, agCs). In addition, the spatial 
links appearing in the fermion action are stout-smeared 
to reduce radiative corrections. 

The anisotropic lattice breaks hypercubic symmetry 
and introduces two new parameters: the bare quark and 
gluon anisotropics, denoted and which must be 
tuned so that the measured anisotropy ^ — as/at takes a 
target value, which must be independent of the physical 
probe used to perform the measurement. In dynamical 
QCD the tuning procedure is quite complicated but is 
now well-understood and is described in ref. 33]. 

The static quark is simulated with the Eichten-Hill ac- 
tion [l3] and the static quark propagator is then the prod- 
uct of unsmeared temporal links. 

F. Simulation details 

The study was performed at a single lattice spacing 
with Nf — 2. Some details of the simulation, for both vol- 
umes, are in Table HIl The lattice spacing is the same on 



Volume 


P 


atrriq 


as 






Configs 


8^ X 80 


1.508 


-0.057 ' 


-0.17fm 


7.43 


8.42 


232 


12^ X 80 


1.508 


-0.057 ' 


-0.17fm 


7.43 


8.42 


245 



TABLE II: A table of simulation parameters for the two vol- 
umes used in this study. 



both volumes and the light and sea quarks are degenerate 
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in the simulation. The determination of the bare quark 
anisotropy, = 7.43 and the bare gluon anisotropy, 
= 8.42 is described in ref. [s^l- These values were 
chosen to obtain a physical anisotropy ^ = 6. Mea- 
surements of the anisotropy from the light pseudoscalar 
dispersion relation on the small and large volumes yield 
values of 5.94(6) and 5.82(5) res pec tively, which are con- 
sistent with the analysis of ref. [33| . For this parameter 
set, the ratio mi^jmp = 0.54. Our study is therefore most 
applicable to the spectrum of Bs mesons. The temporal 
lattice spacing, determined from the spin-averaged IP- 
IS splitting in charmonium, on the same ensemble [s^ . 
was found to be a^p_;^g = 7.04(3)GeV, which implies a 
value for the spatial lattice spacing, of approximately 
0.17 fm. Simulations were performed on 8"^ x 80 and 
12^ X 80 lattices. The use of two spatial volumes allowed 
us to identify possible finite-size effects and probable mul- 
tiparticle excitations in our results. 

The gauge configurations were generated using the 
standard Hybrid Monte Carlo (HMC) algorithm. A more 
detailed description of the particular improvements to 
the HMC algorithm, specifically for anisotropic lattices 
is in reference [33] ■ The generation of the ensemble of 232 
gauge configurations on the small lattice required approx- 
imately 5000 CPU hours while the 245 configurations on 
the larger lattice were generated in approximately 15000 
CPU hours. 

On each volume, we computed the exact contribution 
of the 50 lowest-lying eigenmodes to the light quark prop- 
agator using the ARPACK implementation of the Lanc- 
zos algorithm to determine the eigenvectors. The choice 
iVe„ = 50 was based on earlier analyses carried out on the 
small volume. For similar run parameters, it was found 
that computing the contribution of a few low-lying eigen- 
modes to the fermion propagator exactly did result in a 
significant improvement in signal resolution compared to 
a purely stochastic estimate for the propagator. However, 
the gain in precision achieved by including more eigen- 
modes quickly decreases so that, for our purposes, it was 
not worth considering greater than 50. To evaluate 
the contribution of the higher modes to the light fermion 
propagator on the small volume, time and colour dilu- 
tion was applied to two independent noise sources. On 
the larger volume, a single noise source was used, which 
was also diluted in time and colour indices. 

Once the light fermion propagators had been com- 
puted, the evaluation of the correlation matrices was rel- 
atively straight-forward. Regarding the use of stout-links 
in the interpolating operators, we found that satisfactory 
results were obtained with 30 iterations of the smearing 
algorithm using a weighting factor of p = 0.025. In order 
to construct correlator matrices, 8 levels of Jacobi smear- 
ing were applied to the light quark fields. The quark field 
after j applications of the smearing algorithm is given by 




where q^^^ (x) is the original unsmeared quark field. On 
both lattices, k was chosen to be 0.1 and was fixed at 
12. For the spatially-extended interpolating operators, 
we varied the lengths of the gauge-covariant paths, test- 
ing offsets which extended over a single lattice spacing 
and two lattice spacings in a given direction. We found 
that the single-link operators coupled optimally to the 
ground states in each lattice irrep and these are the op- 
erators which were used in our final measurements. 



III. RESULTS 



In Figures [2] to [6] we plot the effective masses, given by 



m. 



(") 

effecti' 



In 



Cnn{t)ICnn(t^V) 



(14) 



for the correlation matrices in each of the five represen- 
tations considered. These results on the larger (12'^ x 80) 
lattice show the high quality of the data. Convincing 
plateaux, with small statistical errors, are seen for as 
many as five or six states in each irrep. While the physi- 
cal significance of the higher- lying levels cannot be pinned 
down without further study the very precise data give us 
confidence in our determination of the lower-lying energy 
states. 




(13) 



FIG. 2: The efFective masses in the Giu irrep, containing the 
static-light S wave, determined from the optimised correlation 
matrix. The ground state and 5 excitations are shown. 



A. Fitting Procedures 

The vectors, v„, used to determine the optimised cor- 
relation matrices were obtained by solving the eigenvalue 
problem in Eq.(|3]) on time slices to — 2 and — 10, 
although we tested other time slices to check for consis- 
tency. The resulting correlation matrices are then guar- 
anteed to be diagonal only on time slice 2 and time slice 
10. However, for high-quality data we expect the off- 
diagonal entries of the optimised correlation matrices to 
be very small over some time interval. In this interval, 
estimates for the excited state energies can be obtained 
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FIG. 3: The effective masses in the Gig irrep, corresponding 
to the 1^ P wave. The ground state and 5 excitations are 
shown. 



FIG. 5: The effective masses in the irrep, which contains 



the 



and I D-wave states. The ground state and four 



excitations are shown. 





FIG. 4: The effective masses in the Hg irrep, associated with 
the 



3 + 



2 P wave. The ground state and five excitations are 
shown. 



FIG. 6: The effective masses in the G2u irrep, whose lowest 
lying state is a D wave. The ground state and five excitations 
are shown. 



by fitting to the diagonal entries in the correlation ma- 
trices. The effective masses in Figm^es [5] to IHl vifhich are 
defined in terms of the diagonal entries of the correlation 
matrix, demonstrate that this assumption is reasonable. 

However, the most reliable estimates for the low-lying 
energies are found by including off-diagonal matrix en- 
tries in the fit. To implement this, submatrices involving 
operators which project onto a smaller number, M, of 
the lowest-lying states are extracted from the full 8x8 
optimised matrices. The energies of these states En are 
determined by fitting to the entries of this submatrix 
CpQ{t) using the ansatz 

CpQ{t) = Z^zQe-^"\ (15) 

On the small volume, we fit to a 3 x 3 submatrix in the 
Giu irrep, and to 2 x 2 matrices in the other irreps. The 
enhanced signals on the larger volume meant that it was 
possible to fit to a 4 X 4 matrix in the Giu irrep and 
3x3 correlation matrices in the other representations. 
All the fits use a correlated ^^-minimisation algorithm 
and statistical errors are determined from 1000 bootstrap 



samples. 

Since the fitting routine allows for the fact that the op- 
erators v„ may not be orthogonal, it permits us to fit to 
time ranges much greater than the optimisation time to ■ 
The optimal fit ranges were determined from a sliding 
window analysis where t^ax was fixed at an appropriate 
value and imm was varied; the stability of the fitted val- 
ues and the goodness of fit for varying tmin were then 
examined. Sample sliding window plots, for the S wave 
irrep and the Hu irrep, are shown in Figures [7] and [5] 
The fitted masses are seen to be consistent over a range of 
values of tmin with a good x^/Nd.f . Similar results were 
found for the five irreps considered. The final fit ranges 
and the best-fit energies for each of the lattice irreps are 
given in Table IIIII In each of these fits, we disregard 
the result for the highest energy included in the fit since 
we expect that result to suffer from contamination from 
higher levels. As an example, for the S-wave (Gi„) irrep 
we believe the ground and first two excited states are con- 
vincingly determined in our fit procedure. The signal for 
the remaining higher-lying excitations is still remarkably 
good, as shown in the effective mass plots and we esti- 
mate the energies of these states using single and double 
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S 0,25 - 



FIG. 7: A sliding window plot for the Giu irrep. The max- 
imum time slice for these fitted masses is tmax = 30. The 
fitted masses in the ground and first three excited states are 
stable, with good x^/Nd.f., when tmin is varied in fits to the 
full correlation matrix using four exponentials. Although the 
range of stable values of tmin is smaller than seen in Figure |8] 
the number of states determined is larger. 



Irrep {t 


min, ^max. 


Channels 


xVNd.f. 






ground I''* 2"'' 




Giu 


(24,30) 


0.1379(2) 0.210(1) 0.255(3) 


1.15 


Gig 


(25,30) 


0.1906(7) 0.246(4) 0.291(6) 


1.37 


Hg 


(25,30) 


0.1956(7) 0.257(2) 


1.02 


Hu 


(22,35) 


0.239(2) 0.283(3) 


1.00 


G-iu 


(24,30) 


0.230(4) 0.283(8) 


1.35 



TABLE III: The best-fit energies of the ground states and 
higher excitations of each of the five irreps considered. The 
energies are given in units of a^^ ■ The fit ranges used and 
the x'^/Nd.f. values are also given. These results are obtained 
on the larger of our two volumes, 12"^ x 80. The errors are 
statistical only. 



correlation matrices on each time slice, as proposed in 
ref. [3^ and found that both approaches agree over the 
fit intervals. 
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FIG. 8: A sliding window plot for the irrep. Here tmax = 
35 and the ground and first excited states are seen to be stable, 
with good x^/Nd.f., with respect to changes in Wn. 



exponential fits to the relevant correlation matrices. 
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(16) 
(17) 



All the fits were done using a correlated ^^-minimisation 
algorithm. The statistical errors on the fits were again 
estimated from 1000 bootstrap samples. Using the dou- 
ble exponential form it was possible to fit to relatively 
early time ranges, while the single exponential was ap- 
plicable to later time intervals. We have checked that 
fitting these ansatze over different time intervals yields 
consistent results for all the lattice irreps. In addition we 
have applied this fitting procedure to the lower-lying en- 
ergies and found that the resulting fitted energies agree 
within errors with the energies determined using the full 
correlation matrix, as described above. 

We have also compared the effective masses defined in 
Eq. to effective masses obtained by diagonalising the 



B. The spectrum 

In the static limit, hadron energies receive an unphys- 
ical contribution from the static quark self-energy. How- 
ever, this contribution cancels in energy differences which 
are therefore of direct physical significance. The energy 
differences between the ground state in Gi„, which is the 
S wave, and higher states are shown in Figure [HI The 
solid bars are those states determined by a fit to the 
full correlation matrix. In this case we are satisfied that 
these numbers correspond to the eigenstates of the lat- 
tice Hamiltonian. The dashed boxes denote energy levels 
obtained by fitting to the diagonal elements of the cor- 
relation matrix as described earlier. Although the data 
indicate that the dominance of the diagonal entries is 
a reasonable assumption we nevertheless present these 
results as less conclusive. Table IIVI lists the splittings 
shown by the solid bands in Figure [9] in physical units. 
The lattice spacing has been set from the spin-averaged 
IP-IS splitting in charmonium determined on the same 
configuration set and described in ref. [34| . 



C. Finite volume effects 

The first indication that finite volume effects might 
be significant comes from the slight discrepancy between 
the quark anisotropy on the small and large volumes. 
The small lattice used in this study has a spatial vol- 
ume of about (1.35 fm)'^. The spatial volume of the large 
lattice is (2.03 fm)'^. Comparing the lowest energy lev- 
els in each of the irreps across both volumes, we find 
small but significant shifts in the energies of the Giu, 
Gig and Hg states. The lowest energy levels in the 
and G2u representations remain constant within statis- 
tical errors. The Gi„ energy level decreases slightly on 
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H 



H 



Irrep Energies 

8^ X 80 12^ X 80 



- 


Giu 


0.1411(5) 0.1379(2) 




Gig 


0.174(1) 0.1906(7) 




Hg 


0.183(2) 0.1956(7) 






0.233(3) 0.239(2) 




G2u 


0.224(3) 0.230(4) 



TABLE V: A comparison of the fitted energies of the lowest- 
lying states in each irrep on two volumes. These values are 
determined in similar fits to the full correlation matrix. The 
largest effect is seen in the P wave states. The energies of 
these states increase on the larger volume while the splittings 
between them decrease. 



FIG. 9: The splitting between the ground state S wave 
(ground state in the Giu irrep) and the other states deter- 
mined in this study. We show the splitting between states 
determined from fits to the full correlation matrix as solid 
bands. Such fits ensure that even small off-diagonal elements 
of the correlation matrix, which represent leaking between the 
levels, are taken into account. The shaded bands are those 
splittings determined, for the higher-lying states, using fits to 
the diagonal elements of the correlation matrix only. As dis- 
cussed in the text the data indicate that this is a reasonable 
procedure, nevertheless we present these states as less con- 
clusively determined. The possible inversion of ordering seen 
in the D-wave orbital is discussed further in Section |V] of the 
text. 



channel 


AE (GeV) 


Giu, 1st excitation 
Giu, 2nd excitation 


0.504(8) 
0.82(2) 


Gig, ground state 
Gig, 1st excitation 


0.371(6) 
0.76(3) 


Hg, ground state 
Hg, 1st excitation 


0.405(6) 
0.84(2) 


Hu, ground state 
Hu, 1st excitation 


0.706(2) 
1.02(2) 


G2u, ground state 
Giu, 1st excitation 


0.65(3) 
1.02(6) 



TABLE IV: The splitting between the ground state in the 
Giu (S-wave) irrep and the other states determined in this 
analysis by fits to the full optimised correlation matrices. The 
lattice spacing was determined from the spin-averaged IP- IS 
splitting in charmonium, on the same ensemble. 



the larger volume. The energies of the Gig and Hg states 
show the strongest volume dependence. The energies of 
these states increase on the large volume and the splitting 
between them decreases. This is illustrated in Table |V] 
which shows the fitted masses, for the lowest-lying state 
in each irrep on two volumes. These states are expected 
to correspond to P wave excitations of the static-light 
meson. Although a noticeable volume dependence in en- 
ergies can be a signal for scattering states involving par- 



ticles with non-zero spatial momentum, one expects the 
energies of such states to decrease on a larger spatial vol- 
ume. Furthermore, unless the ground state energies in 
Hu and G2u also correspond to multiparticle states, such 
an assignment is highly unlikely. The possibility that the 
observed energy levels correspond to two-particle scatter- 
ing states will be discussed in further detail in Section. IIVI 



D. Quark mass and lattice spacing effects 

The results presented here have been obtained on 
gauge configurations generated with a probability density 
which includes the effects of two degenerate sea quark 
flavours. The mass of the sea quarks was close to the 
strange quark mass and the same value was used for the 
mass of the light valence quark. Therefore, we have sim- 
ulated a well-defined unitary field theory and our results 
are free of quenched pathologies. However, we have not 
explored the inclusion of additional flavours or the use of 
lighter quark masses. Recent developments in algorithms 
indicate that these improvements can be implemented 
with moderate computing resources using the same ac- 
tions in the near future [36l|. 

The largest uncertainty in our results is due to the 
fact that we have carried out simulations at just a sin- 
gle lattice spacing. Therefore, although we are confi- 
dent that the energies we have computed correspond to 
eigenstates of the lattice Hamiltonian, we are unable to 
extrapolate these results to the continuum limit. More- 
over, the spatial lattice spacing used in the simulations 
was quite coarse, Os ~ 0.17 fm, and cutoff effects in our 
results may be large although we hope that the use of 
improved actions keeps these errors under control. 



IV. MULTIPARTICLE STATES 

It is possible that a number of the energy levels we have 
computed are not the energies of bound-state resonances 
but correspond to multiparticle states. The pion mass 
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measured on this configuration set is atruj^ ~ 0.05579 
and one can expect two-particle states consisting of a 
static-light meson and a pion in particular to contribute 
to the energy regime shown in Figure [5] Since we have 
obtained results on two spatial volumes, it should be 
possible to identify multiparticle resonances by compar- 
ing the elements of the optimisation vectors v„, which 
are computed in the variational method, on both vol- 
umes. The normalisation condition for these vectors is 
v5^C(to)v„ = 1. A vector projecting onto a bound state 
is, to a good approximation, volume independent. On the 
other hand, a vector which isolates a multiparticle state 
shows a strong, well-defined volume dependence [37| . 
This dependence comes from the fact that, for example, 
a state consisting of two weakly-interacting particles in 
a finite spatial volume has a spectral weight which scales 
with the inverse volume. The ratio of the volumes used 
in this study was 3.375, implying that multiparticle ex- 
citations should be readily identifiable. However, studies 
of possible pentaquark candidates have shown that con- 
siderable care must be taken when trying to separate 
single-particle energies from the energies of scattering 
states on the basis of the volume dependence of spec- 
tral weights [H, [H, [13] . In particular, when attempting 
to isolate the spectral weight of a single state one needs 
to check for contamination from higher states. We there- 
fore consider the volume dependence of the eigenvectors 
in conjunction with fit values for the overlap of the opti- 
mised interpolating operator onto the state under consid- 
eration. We examine the eigenvectors which project onto 
the three lowest states in the Giu irrep and the two low- 
est states in each of the other irreps, for which we have 
reliable energy measurements. However, we find no evi- 
dence for multiparticle contributions to the energy spec- 
trum shown in Fig. [9l although such states undoubtedly 
exist. We therefore conclude that the interpolating op- 
erators used have negligible overlap with low-lying mul- 
tiparticle resonances. Perhaps this is not so surprising, 
by construction one expects the interpolating operators 
to couple strongly to just a few low-lying single-particle 
states. In fact, given that the operators used in the vari- 
ational bases differed only in the smearing of the quark 
fields, it is rather remarkable that the variational method 
appears to yield good results for so many single-particle 
resonances. 

Ultimately, the issue of multiparticle states in the 
static-light spectrum will be resolved by using much 
larger operator bases incorporating interpolating oper- 
ators specifically chosen to project onto multiparticle 
states. However, using the data available to us it is possi- 
ble to make some simple predictions for the multiparticle 
spectrum. If we assume that the interaction between par- 
ticles in a scattering state is very weak, we can determine 
approximate values for the threshold energies for multi- 
particle states in each of the irreps. For a large enough 
spatial volume the lowest-energy two-particle states con- 
sist of the ground-state S-wave static-light meson and a 
pion which may have non-zero momentum. We use peri- 



odic boundary conditions in the spatial directions so the 
components of the pion momentum p are quantised in 
units of 2TT/asL. A pion at rest transforms according to 
the Aiu irrep of Oh- However, by applying the elements 
of Oh to a pion field with non-zero momentum, one can 
generate more complicated single- valued representations 
of the group. Table IVII gives the decomposition of pion 
representations into their constituent irreps for the small- 
est values of |p| allowed by the lattice. Representations 



p 


Irreducible content 


(0,0,0) 




(1,0,0) 


Aiu ®Eu® Tig 


(1,1,0) 


e K e Tig e Tag e T2„ 



TABLE VI: An example of the decomposition of pion repre- 
sentations of Oh, labeled by the pion momentum in units of 
2-K jasL, into their constituent irreps. 

for the two-particle scattering states of interest are given 
by the direct product of the pion representations with 
G\u- We then decompose these representations in terms 
of the double- valued irreps of Oh- By considering differ- 
ent pion momenta one can deduce the two-particle scat- 
tering content of each of the irreps. This procedure is 
in complete analogy to the determination of the angular 
momentum content of the irreps described in section lTl Al 
The energy of a scattering state is approximated by the 
sum of the energies of the decoupled states. 

S2P » Es.L. + ^m2+p^ (18) 

where mi is the rest mass of the light meson appearing 
in the scattering state. The threshold for multiparticle 
states in a particular irrep is the minimum such energy 
appearing in that representation. 

The analysis is complicated by the fact that on the 
12'^ X 80 lattice the smallest non-zero momentum compo- 
nent is approximately 0.088, in units of a^^, which lies 
between the pion and rho masses. Therefore, the lowest- 
energy scattering state in a given irrep may contain a rho 
meson at rest rather than a pion with non-zero momen- 
tum. Moreover, the unit of momentum is also of the same 
order of magnitude as the splittings shown in Figure[9l so 
it is possible that the threshold levels in some irreps cor- 
respond to states containing a static-light meson which 
lies in a representation other than Gi„. We have con- 
sidered both possibilities and find that on this spatial 
volume the threshold energies do in fact correspond to 
two-particle scattering states involving only the S-wave 
static-light meson and a pion. However, the threshold 
energy in Hg is degenerate within errors with the energy 
of a scattering state consisting of an S-wave static-light 
meson and a rho meson at rest. The lowest threshold 
across the irreps is in G\g and it is associated with a 
pion at rest. The threshold energies in Giu, Hg and 
correspond to a pion with momentum (1,0,0) in units 
of 2t: /asL- In the lowest allowed pion momentum 
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is (1,1,0). In Figure [TPl we have inserted the thresholds 
for muhiparticle excitations into the sphtting plot previ- 
ously shown in Figure [HI The dashed red lines indicate 



FIG. 10: Splittings between the energy of the ground-state 
S-wave static-light meson and other low-lying measured ener- 
gies. The dashed red lines show the scattering thresholds in 
each of the irreps. The position of the threshold in the Gig 
irrep corresponds to the rest mass of the pion. The thresholds 
in the other irreps correspond to two-particle states contain- 
ing a pion with non-zero momentum. 

the scattering thresholds, and we see that multiparticle 
states should be detectable in all five irreps over the range 
of energies investigated in this study. 



V. DISCUSSION 

Although a final determination of the excited state 
spectrum requires simulations at several lattice spacings, 
it is possible to draw some conclusions from the results 
obtained here. 

First we consider energy levels in Gi„. The lowest- 
lying state in this irrep is the Jp = | S wave. The 
next lowest angular momentum value allowed in this ir- 
rep is = |. The I state corresponds to an L=4 or 
G wave excitation. We expect such a state to be ex- 
tremely heavy and also to appear as a degenerate level 
in other lattice irreps. Since this is not the case we there- 
fore conclude that the energy-levels lying directly above 
the ground state correspond to radial excitations of the 

i . Similarly, the lowest-energy state in the Gig irrep 
1 + 

is the ^ , and we expect the first excited state to be a 
radial excitation with the same quantum numbers. The 
interpretation of the higher-lying states is slightly more 

7 + 

ambiguous due to the appearance of the ^ F-wave state 
in this irrep. In the Hg irrep the ground state has quan- 
tum numbers Jp = This state lies slightly above 

the ground state in the Gig irrep. These are P-wave 
excitations. The data show solid evidence that the natu- 
ral ordering prevails in the P-wave sector (m3_|_ > mi^). 
While this splitting is very small it has been measured in 
this work to be different from zero at the 5(t confidence 



level. The discussion of natural ordering is made more 
complicated by the fact that the Gig ground state lies 
close to the estimated S wave decay threshold. The first 
excited state in Hg lies just above its counterpart in Gig 
and we also identify it as a radial excitation with P wave 
quantum numbers. 

To make sense of the energy levels in i7„ one needs to 
consider them together with the states appearing in G2u- 
To begin with, we will assume that the variational analy- 
ses have succeeded in determining all the low-lying single 
particle energies in the ff„ and G2u irreps, and consider 
how the relative ordering of D-wave states should effect 



the simulation results. Recall that both the | 



and I 



appear in the Hu irrep and only the | D-wave con- 
tributes to the G2u irrep. In the standard ordering of 
states, the | is lighter than the | which implies that 
the lowest-lying state in Hu should lie below the lowest- 
energy state in the G2u irrep. One of the higher energy 
levels in the H^ irrep, corresponding to the | multiplet, 
should then be degenerate with the lowest energy level in 
G2u- This scenario is clearly not supported by the data. 

Another possibility is that inversion of the D-wave mul- 
tiplets does occur. In this case, the ground state energies 
in Hu and G2u are degenerate up to lattice artifacts. The 
energies of the first excited states in these irreps overlap, 
and these states may correspond to a radial excitation of 
the I . This would mean that the | state is the third 
level in the Hu irrep which has no counterpart in the G2u 
irrep. However, it seems very unlikely that the splitting 
between the D-wave multiplets is so large. It could also 
be the case that systematic errors in the data are partic- 
ularly severe in these channels, and that at finer lattice 
spacings and larger volumes the second energy level in 
Hu shifts towards the ground state while the third en- 
ergy level aligns itself with the second level in the G2u 
irrep. One could then identify the second level in Hu with 
the I excitation, and the third level in 7J„ and the sec- 
ond state in G2u would correspond to a radial excitation 
of the I . However, given that the results were obtained 
using improved actions on a (2.03 fm)'^ spatial volume, 
it is improbable that systematic errors could account for 
such significant shifts in the spectrum. 

The inconsistencies between the data and the possi- 
ble orderings of states seem to suggest that the assump- 
tion on which the preceding argument is based is incor- 
rect, and we have not succeeded in measuring all of the 
lowest-lying single-particle energies in Hu and G2u- In 
particular, there appear to be gaps in the spectrum of 
the Hu irrep. This is not surprising when we recall that 
each variational basis consists of operators which have 
identical offsets. It appears that a complete determina- 
tion of the spectrum will require a much greater variety 
of basis operators with varying degrees of overlap with 
all the low-lying states. Naively, the fact that the two 
lowest-lying measured energies in Hu are close to their 
counterparts in G2u might indicate that these are the 
energies of | states. However, we have already noted 
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that the third level in Hu does not have a partner in G2u, 
and there is no reason why the lower energies in _ff„ could 
not correspond to | excitations. Therefore, identifica- 
tion of the lowest measured energy levels in the Hy, irrep 
requires further investigation. 

In spite of the difficulty in interpreting the data, our 
study of the and G2u irreps has been quite successful, 
and we would like to emphasise two points. Firstly, we 
have been able to measure the energies of the ground- 
state I D-wave and a radial excitation. These corre- 
spond to the first and second energy levels in the G2u 
irrep. Secondly, the possibility that the interpolating op- 
erators are not coupling to all the low-energy states has 
only come to light because we have been able to precisely 
determine a number of energies in each of the irreps which 
we have compared to theoretical predictions. 

Finally, it is worth comparing our data with experi- 
mental data for the Bg spectrum. The energy difference 
between the _B*j(5850), with unknown quantum num- 
bers, and the ground-state pseudoscalar is approxi- 
mately 485 MeV. This number is very close to our value 
for the splitting between the ground-state S wave and its 
first radial excitation. However, once l/mq corrections 
are taken into account, we see that the splitting is also 
consistent with a P-wave identification for the heavier 
state. 



spatially-extended interpolating operators which project 

onto the excited states. The construction of these op- 
erators and proper interpretation of the final numerical 
results depend on a clear understanding of the spatial 
symmetries of the lattice, and our results are given in 
terms of the irreducible representations of the octahedral 
point group. 

We have succeeded in determining a number of excited 
state splittings in five of the six double- valued representa- 
tions. A comparison of results obtained on two different 
spatial volumes indicates that that our results correspond 
to the energies of single- particle bound states. We have 
presented solid evidence (at the 5(t level) that the natu- 
ral ordering of P-wave static-light mesons prevails. We 
have been able to identify a number of radially excited 
states however a complete survey of the low-energy spec- 
trum will require further work including, ultimately, an 
extrapolation to the continuum limit. 

A natural progression in this study would include the 
use of interpolating operators with more complicated off- 
sets which might better couple to heavier single-particle 
excitations, and the use of operators specifically con- 
structed to couple to multiparticle states. We are also 
confident that a similar approach to the one used here 
can shed considerable light on the spectrum of heavy- 
light baryons. 



VI. CONCLUSIONS 

We have presented an ab initio study of the spectrum 
of static-light mesons in TV/ = 2 lattice QCD. Previously, 
computational constraints have made it difficult to deter- 
mine the excited state spectrum to satisfactory accuracy. 
We have overcome these constraints using an efficient es- 
timate for all elements of the lattice Dirac propagator. 
With this method and exploiting translational invari- 
ance, we are able to average correlation functions over the 
whole volume of the lattice, yielding a dramatic reduction 
in the variance of our Monte Carlo estimates. The use of 
all-to-all propagators also facilitates the construction of 
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